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ABSTRACT 

Aims. What is the origin of the large-amplitude variability in Wolf-Rayet WN8 stars in general and WR123 in particular? A dedicated 
spectroscopic campaign targets the ten-hour period previously found in the high-precision photometric data obtained by the MOST 
satellite. 

Methods. In June-August 2003 we obtained a series of high signal-to-noise, mid-resolution spectra from several sites in the AA 4000 - 
6940 A domain. We also followed the star with occasional broadband (Johnson V) photometry. The acquired spectroscopy allowed a 
detailed study of spectral variability on timescales from ~ 5 minutes to months. 

Results. We find that all observed spectral lines of a given chemical element tend to show similar variations and that there is a good 
correlation between the lines of different elements, without any significant time delays, save the strong absorption components of 
the Hei lines, which tend to vary differently from the emission parts. We find a single sustained periodicity, P ~ 9.8 h, which is 
likely related to the relatively stable pulsations found in MOST photometry obtained one year later. In addition, seemingly stochastic, 
large-amplitude variations are also seen in all spectral lines on timescales of several hours to several days. 

Key words. Stars: individual: HD 177230 = WR 123 - Stars: Wolf-Rayet - Stars: winds, outflows - Stars: oscillations 



^ , 1. Introduction runaway frequency is relatively high, and in some cases their 

k> _, . . TTI7 ,„ m . r , „„„ , , , distance from the Galactic plane is la rge (Moffat & Isserstedtl 

^ . Population I Wolf-Rayet (WR) stars of the WN8 subclass have 1980; Moffat 1989; Moffat gt aJ 199g) ^ ^ ^ hoto _ 

S ' b6en * 6 u SUbjeCt ° f ™ ous , studle « 0Wln g to their fug- metriCj po i arimetric and spectroscopic monitoring campaigns of 

matic behavior. These direct descendants of massive and hot WN8 gtars haye consistently shown the highest level of intrin- 

Of stars are markedly distinct from other WR stars. Indeed, sic variabmty among all WR stars ( . Moffat & shara 1986; 



no member of this class is known to have a close massive O- Lamontagne & Moffat 1987; Drissen et'al. 1987; Balona et aL 



type companion. They avoid star clusters and associations. Their 19g9; Robertetal . 1989; R ' obert 1992; Antokhin et al. 1995^ 



* Based on observations obtained at the Canada-France-Hawaii |Marchenko et al.||1998|). However, despite the claims of multi- 
Telescope (CFHT) which is operated by the National Research Concil pie photometric periods ranging from hours to days, there is no 
of Canada, the Institut National des Sciences de l'Univers of the Centre period that could be labeled as unique to all the campaigns for 
National de Recherche Scientifique of France, and the University of any given WN8 star. Discovery of any stable periodicity pattern 
Hawaii. Based also on observations made with ESO Telescopes at the in a WN8 star would help to establish whether the star is a binary 
La Silla Observatory, under programme ID 271.D-5025. or a s i ng i e rotator with an inhomogeneous wind or, alternatively, 
** Photometric date presented in Figure 1 are avail- to refine the rather controversial pulsation models of WR stars, 
able in electronic form at the CDS via anonymous ftp 

to cdsarc.u-strasbg.fr (138.79.128.5) or via All of the above applies specifically to our target, WR 123 

http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A+A/ (see Marchenko et al. 1998 for references), which shows an 
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extremely complex, rapidly evolvi ng light curve, along with 
seemingly aperiodic polarimetric dSt-Louis et al.l Il988l) and 
high-amplitude line-pro f ile variability (Masse y & Contil 11980; 
lLamontagne et al.l 119831: lMarchenkoetal.lll998l) . All searches 
for periodic signals provided rather inconsistent results until 
2004, when an inte nse, long-term, broa d-band photometry cam- 
paign of WR 123 dLefevre et al.1 12005) was performed by the 
MOST (Microvariability and Oscillations of STars) satellite. It 
was shown that large-scale, short-lived (typically, less than 1 
week) photometric fluctuations were dominated by the ever- 
variable continuum. Only one relatively stable, small-amplitude 
period was present in the entire data set at P « 9.8 h. Does this 
period indicate the detection of a second spectroscopic binary 
system among WN8 star£] or does it provide a reliable observa- 
tion of periodic stellar activity in a WR wind, like pulsations or 
co-rotating density structures? 

While broadband photometry documents the variability of 
the (mainly) stellar continuum emanating from the inner WR 
wind, spectroscopy probes the layers outside it. Moreover, since 
the emission lines of different chemical species and/or different 
ionization potentials are generally formed in different regions of 
the wind, the timing and correlation of profile variations detected 
in such lines may help for determining how perturbations prop- 
agate in a stratified wind. In this paper we present the results of 
intense spectroscopic monitoring of WR 123. 

2. Observations and data reduction 

Originally, we planned to obtain spectroscopic monitoring in 
parallel with the MOST run, initially scheduled for June/July 
2003. Therefore, we carried out spectroscopic observations of 
WR 123 during 5 runs on 4 different telescopes (Table 1). 
Unfortunately, the launch of MOST was delayed by one year. 
A small-scale spectroscopic campaign run in parallel with the 
2004 MOST observations proved to be of critical importance 
and provided a picture, though rather sketchy, of line-profile 
variability (lpv) supe rposed on an ever-changing continuum flux 
dLefevre et alJ l2005). Here, working with a unique, much more 
extensive data-set, we intend to obtain a more detailed descrip- 
tion of the lpv, albeit outside the MOST campaign. The acquired 
spectroscopy allowed us to follow lpvs on timescales from ~ 5 
minutes (NTT spectra only) to months (all spectra). 

In parallel with the 2003 spectroscopic monitoring, we 
also performed broadband (Johnson V) photometry using the 
24-inch telescope of the Bell Observatory (Western Kentucky 
University, USA). The differential photometry was carried out 
using the weighted combination of comparison stars in the 4.4 x 
6.6 arcmin CCD field. We used comparisons with the following 
(epoch 2000) RA and DEC coordinates: (a = 19:03:58.2; 6 = - 
04:18:48),(a = 19:04:04.6; 6 = -04:17:35), (a = 19:03:45.3; 
6 = -04:19:30), (a = 19:03:54.7; 6 = -04:20:48), and 
1.0:0.75:0.25:0.25 weights, as well as a star at (a = 19:04:13.2; 
6 = -04:20:3 1) as a control star. We plot the available photomet- 
ric data and corresponding + 1 cr errors in Figure 1 . All the spec- 
troscopic and photometric data were processed using standard 
irai0 software packages. Special care was taken for the normal- 
ization of the spectra. First, a mean was made for each run. Then 



1 The runaway WR148 = HD197406 is the first detected WN8 
spectroscopic binary with a moderate-mass, unseen companion 
dMarchenko et al.f l996) 

2 iraf is distributed by the National Optical Astronomy Observatory, 
which is operated by the Association of Universities for Research in 
Astronomy, Inc. under cooperative agreement with the National Science 
Foundation. 



each spectrum of a run was divided by the run mean and the ratio 
fitted with a low order Legendre polynomial (between 4 th and 8 th 
order). The original individual spectrum was divided by this fit, 
and was therefore at the same level as the run mean. When this 
procedure was done for each run, the run means were then put at 
the same level by using the same procedure as described above. 
This allowed us to put all individual spectra at the same level. 
Then, we combined all run means into a global, high quality 
mean, which was then fitted in selected pseudo-continuum re- 
gions, i.e. wavelength regions where large emission-lines do not 
dominate. These regions are listed in Table 2. To avoid having 
to use high order polynomials for the fit of a very large wave- 
length region, we cut the spectra in 5 different sections to fit, i.e. 
3999 _ 4454 A, 4400 - 4794 A, 4750 - 5294 A, 5229 - 6032 A 
and 5999 - 7000 A. Finally, the fitted continuum function is ap- 
plied to each individual spectrum. The error on the continuum 
normalization measured as the standard deviation of individual 
spectra around the continuum function is typically of 0.5 %. 
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Fig. 1. The contemporaneous broad-band (Johnson V) photom- 
etry of WR 123 obtained at the Bell Observatory. AV is the 
difference between the magnitudes of WR 123 and comparison 
stars (see text). Each measurement is plotted with corresponding 
± 1 cr error bars. The upper panel shows all acquired data while 
the lower panel depicts two nights from August 2003. 



3. Results 

3. 1 . Temporal variance spectrum 

In order to identify the variable parts of the spectrum of WR 1 23, 
we adopt th e "temporal variance sp ectrum" (TVS) approach de- 
veloped by Fullert on et al.l (119961) . To minimize the influence 
of intrinsic radial velocity (RV) variations of WR 123 on the 
TVS, we co-align the spectra using cross-correlation. All the 
spectra were resampled after the co-alignment at the resolu- 
tion of the spectra obtained at NTT, since they are the most 
numerous in our dataset and they have the lowest spectral res- 
olution. The resampling was made using a sine function. This 
method tends to slightly broaden the line profiles (negligibly in 
the case of WR spectral lines), but it does not modify the level of 
noise in the spectra and, therefore, does not affect the TVS. We 
find that all the spectral lines in the observed spectral domain, 
AA 4000 - 6940 A, show lpvs significantly above the adopted 
99 % threshold. All reasonably isolated spectral lines produce 
a similar shape in the normalized TVS plot. However, in a given 
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Table 1. Spectroscopic Observing Campaigns 





La Silla 


Observatoire de 


Canada-France- 


Observatoire du 




Observatory 


Haute-Provence 


Hawaii Telescope 


mont Megantic 


HJD (- 2452800) 


15.7,40.8 -45.8 


28.5 -34.5 


60.9 -64.8 


73.6 


Telescope 


3.58m (NTT) 


1.52m 


3.6m 


1.6m 


Instrument 


EMMI 


Aurelie 


OSIS 


OMM-Spectro 


No. of spectra 


290 


66 


17 


2 


Resolution, A 


5.1 (6 pix.) 


1.4 (3 pix.) 


1.5 (3 pix.) 


1.6 (3 pix.) 


Spectral range, A 


4000-6940 


5225-6100 


4400-6100 


4500-6100 


Signal-to-noise 


>200 


100 


200 


100 



Table 2. Wavelength intervals of the selected pseudo-continuum regions 



wav. int. wav. int. wav. int. wav. int. wav. int. 



(A) (A) (A) (A) (A) 



4001 


-4006 


4495 


-4498 


5095 


-5105 


5620- 


5635 


6197 


-6201 


4044 


-4047 


4585 


-4589 


5245 


-5255 


5760- 


5766 


6335 


-6347 


4268 


-4279 


4740 


-4755 


5372 


-5381 


5990- 


5998 


6638 


-6643 


4420 


-4430 


4830 


-4840 


5515 


-5525 


6021 - 


6025 


6771 


-6794 


4422 


-4428 


4965 


-4975 


5555 


-5565 


6057- 


6061 


6927 


-6933 



P Cygni profile the absorption displays a single narrow peak in 
the normalized TVS spectrum, and the emission part shows a 
wider and more structured peak with stronger variability in the 
central parts (e.g., Figure 2 with Hen A 5412 A and Hei A 5876 A 
as typical examples). This indicates that the detected lpv are not 
directly related to continuum variability. The latter would create 
a predictable global pattern with a shape following the emission- 
line contours. 

3.2. Line-profile and equivalent-width variability 

The TVS analysis shows that, in absolute terms, the most vari- 
able and [relatively] isolated lines are Hei A 4471, 4922, 5016, 
5876, 6678 A; Hen A 4686, 4860, 5412, 6560 A; and NmAA 
4640,41 A. 

To study the nightly changes in the line profiles, we single 
out two strong and least blended transitions, Hen A 5412 A and 
Hei A 5876 A. The most variable part of the Hei A 5876 A line 
is its absorption component, though the emission also changes 
from night to night. No apparent recurrence or coherent behav- 
ior in the variability patterns can be found at first glance. The 
Henvl5412A line does not precisely follow the changes in Hei A 
5876 A. However, the lines behave in a roughly similar manner, 
i.e. the amplitudes of the absorption and emission components 
vary strongly and seemingly randomly on a timescale of days 
and longer. 

From inspection of the remaining lines we conclude that 
what determines the global variability patterns of a line on a 
timescale of days is the strength of the P Cygni absorption com- 
ponent. Indeed, regardless of their species, all lines with a rel- 
atively weak absorption component change similarly to Hen A 
5412 A while all lines with a strong P Cygni absorption behave 
like Hei ^5876 A. 

On a timescale of hours all the lines show similar lpv in 
shape and in relative intensity, if one accounts for spectral blend- 
ing effects. Within a given night, intensity and width can vary 
from 5 % to 30 %, with the latter, if applied to strong emission 
features, in clear excess of the variability expected from the ever- 
changing continuum level on the same timescale. 

In an attempt to quantify the changes for each of our se- 
lected lines, we examine their equivalent widths (EW). Since 



the variability of the absorption and emission components of 
P Cygni profiles do not appear to be correlated, we have treated 
their EW measurements separately by fitting two gaussian func- 
tions, one with a negative intensity and displaced from the 
line center (absorption) and another with a positive intensity 
(emission). The measu r emen t errors were obtained following 
Vollmann & Eversber i d2006). Though realizing that such an 
approach may not be applicable to all considered cases, we note 
that most of the 2-component fits have a reasonably low y 2 an d 
produce fairly good results for the blue-side portions of profiles. 
In an attempt to sample a significant part of the WR wind we 
choose the Hei A 4471 A (which gives better fits than Hei A 
5876 A) and Hen A 5412 A lines, both of which are strong, iso- 
lated and relatively easy to fit. We find that, in general, a stronger 
emission corresponds to a stronger absorption, and vice versa. 
However, for the Hen A 5412 A line, there are frequent dispari- 
ties between the variability patterns of the emission and absorp- 
tion components. Moreover, even if both lines (Hei and Hen) 
vary with comparable amplitudes on a timescale of days, their 
variability within one night can be completely different. Taking 
into consideration the rest of the measured lines, we find that the 
EW of the emission parts of the Hen, Nm and Niv profiles are 
well-correlated, while quite frequently the Hei lines lack corre- 
lation with any other species and among themselves. This could 
be partially blamed on the presence of strong absorption compo- 
nents with a pronounced tendency to show different lpv patterns. 
The lack of correlation is especially prominent in the long-term 
(days-months). However, a tight Hei - Hen correlation may reap- 
pear in some of the short-term (nightly) data segments. 

3.3. Radial velocity and skewness variability 

Previous RV studies of WR 123 have consistently shown the 
presence of small-amplitude vel ocity variations with K < 
25 km s _1 dMarchenko et al.l [T998). Unfortunately, the NTT and 
CFHT data produced large systematic errors related to specific 
instrumental setups. This leaves OHP as the only reliable source 
of RV. We aligned the OHP spectra via cross-correlation with 
the highest S/N spectrum (the first spectrum of the campaign in 
this case), thus producing a mean OHP spectrum. Then we used 
this mean to perform a cross-correlation with individual spectra 
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5400 5500 5600 5700 5800 5900 6000 

wavelength (A) 



Fig. 2. Mean rectified spectrum (solid line) of WR 123 calculated using all the spectra from all the observing runs between 
53 50 A and 6000 A sp ectral range (common to all runs). The TVS spectrum (dash-dotted line) obtained using the formalism 
of lFullerton et alj (1 19961) is overplotted. The dashed line is the threshold of significant variability at a level of 99 % for the TVS. 
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Fig. 3. The de-trended EW measurements of the emission component (left panel) and of the absorption component (right panel) of 
the Hen A 5412 A line folded in phase with P = 0?400 and an arbitrary phase. Five half-nights of the second NTT run are shown 
along with the corresponding 1 <x error bars. 



one more time and measure their RV. RVs were checked with 
the bisectors of the strongest emission lines. While the bisectors 
were giving bigger errors, they match the results coming from 
cross-correlation and show that all the lines are shifting in time 
with the same fashion. 

On a weekly time-scale, the measured RVs vary by up to 
+ 20 km s _1 , well in excess of the typical errors, cr ~ 2 - 
5 km s _1 (obtained from the fit of the cross-correlation pro- 
file), by gradually decreasing from ~ km s _1 to ~ -40 km s _1 
(with respect to the first spectrum of the campaign). A closer 
look shows that the RVs vary from night to night and can change 
by as much as 10 km s _1 within a ~ 4-hour interval. Are these 
RV variations caused by the presence of a companion or simply 
by large-scale lpv? One possible way to verify this is to examine 
the skewness of the main spectral lines, i.e. Hen A 5412 A and 
Hei A 5876 A. Since these two lines clearly dominate the spec- 
trum (and therefore the flux in the cross-correlation profile) in 
the chosen wavelength range, they are also the main contributors 
to the RV measurements. We find that the skewness varies on a 
short timescale of ~ hours. Hence, the variability of the RVs may 
be induced by large-scale, short-term profile variability, save the 



long-term downslope in the RV curve. It is not possible to re- 
late any particular change in the overall line shape to a particular 
change in RV. With some caution, we conclude that the only 
"genuine", i.e., caused by a possible spatial (orbital?) motion of 
the star, RV variability may occur on a fairly long timescale of 
days or weeks. 



3.4. Photometric variations 

Our light curve (Figure 1) has fairly simi l ar ch aracteristics to 
the light curve published bv ILefevre et al.l (2005). Both have a 
variation amplitude of ~ 15-20 % on the long term and ~ 5 % on 
the short term. Hence, we can assume that nothing major in the 
star and the star's wind has changed between the 2003 spec- 
troscopic campaign (this study) and the 2004 MOST campaign 
(ILefevre et al.ll2005l) and that the two datasets can be compared. 
However, the limited photometric accuracy and the poor time 
sampling of our data compared with the unique MOST obser- 
vations does not allow us to say anything more about WR 123's 
photometric variability. 
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4. Period search 

We performed a period search in the permittable range from 10 
min up to one month, using two complement ary approaches. The 
phase-dispersion minimization algorithm (Stelling werll Il97§h 
(PDM) is well-suited for cases with only a small number of ob- 
servations and highly non-s inusoidal underlying signal. We als o 
used the clean algorithm (IScargldll982L Rober ts et al . 1987), 
which has the advantage of taking into account the unevenness 
of the data sampling. 

We searched for periodic signals in the RV data and found 
that the OHP subset reveals no significant periodicities. Period 
search on the nightly mean EWs of the Hen A 5412 A line pro- 
vides inconclusive results, mainly due to the scarcity of these 
particular (1 point/night) data. Searching for short periods, we 
concentrated on our numerous NTT spectra, first carefully re- 
moving any long-term (days-weeks) trends. We find a period 
of 0?400 + C*?004 (v = 2.50 + 0.06 d" 1 ), which is compati- 
ble with th e range of frequencie s between 2.40 d -1 and 2.50 d _I 
claimed by ILefevre et al.l d2005l) from MOST data , is present in 
the de-trended EW measurements of both Hei (peak at 2.46 d _1 ) 
and Hen lines (peak at 2.5 Id -1 ) and is clearly detected by PDM, 
but only marginally present in the CLEANed PS. The de trended 
EW values of Hen A 5412 A folded in phase with this period are 
shown in Figure 3. 

Since the amplitude of the photometric variations on the 
timescale of a few hours is relatively small compared with the 
accuracy of our photometric dataset and due to our poor time 
sampling, no clear period could be found in the light curve. 



5. Discussion 

Most of the lines in the WR 123 spectrum have a P Cygni profile, 
i.e. they contain a blue-shifted absorption component formed in 
the column of wind between the stellar core and the observer, 
as well as an emission component formed in a (more or less) 
symmetric shell surrounding the core. Since the two components 
produce different variability patterns, we treated them separately 
wherever feasible. 

All the emission components vary with similar shapes and 
relative amplitudes, up to 30 % (peak-to-valley) of an average 
line intensity on a short timescale (hours), and up to 55 % on a 
long timescale (days). There is no delay between the variability 
patterns of the different lines; i.e., the variability seems to occur 
simultaneously in the different line-forming regions of the wind. 

The absorptions have a seemingly higher level of variability, 
up to 50 % of their depth on a short timescale, and as high as 
90 % on a long timescale. However, this difference may come 
from the intrinsic difference between the formation process and 
the different range of intensity of the absorption (1 to relative 
intensity) and the emission (1 to I relative intensity, where I is 
the intensity at the peak) parts of P Cygni profiles. Moreover, the 
variations of the two parts of any given P Cygni profile have sim- 
ilar shapes, indicating that they can be caused by the same phe- 
nomenon. All absorption components vary in concert, without 
detectable delay between different line species. Considering that 
the star shows a consistent level of photometric variability that 
practi cally never excee ds ~ 0.15 mag long -term (days-month: 
e.g. iMarchenko etail 119981; ILefevre et all l2005l our study) 
and rarely exceeds ~ 0.05 mag during any particular night of 
ground-based observations (e.g., Figure 1), such a level of lpv 
cannot be a consequence of the continuum variations alone. We 
note that for WR 123 the continuum flux is the dominating (usu- 



ally up to 80 %) component in the broadband photometric optical 
flux. 

In the following we concentrate on the $4 period found in 
the de-trended EW data. We plot the period-folded lpv of the 
strongest and least blended lines He n A 5412 A, Hei A 5876 A 
and Hen A 6560 A in Figure 4. In the bottom parts of the panels 
we plot average profiles derived from four different phase inter- 
vals, = 0.0-0.2,0 = 0.3-0.5,0 = 0.5-0.7,0 = 0.8-1.0. 
It is apparent that practically the whole wind responds to the pe- 
riodic perturbation coming from the core. In general, the absorp- 
tion components are more involved in the periodic changes, as 
they come from a more localized region of the wind compared to 
the emissions. This does not hold true for Hei/i5876A which has 
the largest line-forming region among the considered lines, thus 
the best ability to "smooth" out the phase-dependent variability 
patterns. 
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Fig. 4. Upper: Lpvs observed in the Hen A 5412 A (left), HeiA 
5876 A (center) and Hen A 6560 A (right) folded with the 0?4- 
period. Long-term trends have been removed from the spectra 
and the remainders were enhanced by a factor of 2 for clarity. 
Lower. Binned profiles over the four phase intervals = 0.0 - 
0.2, = 0.3 - 0.5, = 0.5 - 0.7, = 0.8 - 1.0. 

ILefevre et al.l d2005l) already dismissed the possibility that 
stellar rotation could be at the origin of the 9.8 hour period, 
since this would lead to a rotational velocity greater than the 
breakup speed at the surface of a typical WN8 star (when tak- 
ing /?* ~ 15 R e according to Crowth eret al.ll 19951) . Also, they 
showed that such a short period could not be produced by the 
presence of a (compact) companion revolving around the star, 
since the deduced orbit would be smaller than the stellar hydro- 
static radius of WR 123. Finally, they placed a firm upper limit, 
A < 0.2 mmag for P $ 1 h our, on the am plitudes of strange- 
mode pulsations predicted by Glatzel et al. ( 1999). 

Since the publication of the MOST light curve, two 
independent models involving non-radial pulsations were 
proposed to explain the 9.8 hour period. The first, by 
iTownsend & MacDonaldl d2006l) . examines the stability of I = 
1 and ( — 2 g-modes and finds that the ^--mechanism asso- 
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ciated with a hot Fe-opacity bump within the star can be re- 
sponsible for periods ra nging between 3 an d 12 hours. The sec- 
ond mode l proposed by iDorfi et ahl (120061) contradicts the con- 
clusion of Lef evre et al.l d2005l) concerning strange-mode pulsa- 
tions. Indeed, using an adequate stellar radius for WR 123, they 
showed that strange-modes with a period as high as ~ 10 hours 
can be excited. However, both models exploited stellar parame- 
ters b eyond the usually adopted range: iTownsend & MacDonaldl 
(2006)' used a s tellar radius that is far too small for a WN8 star, 
and Dorfi et al 1 (l2006h used a hydr ogen ab undance which is too 
high for WR 123. More recently, iGlatzell (120081) followed the 
evolution of strange mode instabilities into the non-linear regime 
and showed that they generate consecutive shock waves which 
inflate the stellar envelope considerably and thus increase the 
pulsation periods, giving final periods that are consistent with 
that observed in WR 123. It remains unclear how these pulsa- 
tions will affect the stellar wind (thus the wind-formed spectral 
lines) at different radii. Our data show that the 9.8-hour spectral 
variability involves all prominent lines in the observed spectral 
region, thus globally affecting the wind. 

Finally, on a speculativ e note inspired by the suggestion 
of iFoellm i & Moffat (2002), one might wonder if a Thorne- 
Zytkow-like scenario could be operating in WR123 and other 
WN8 stars a nd exciting stellar puls ations. Thorne-Zytkow 
Objects (TZO, iThorne & Zvtkowll 19771) are assumed to be red 
supergiants with a degenerate neutron core. They are formed via 
merger of a massive star and a remnant of a supernova explosion 
in a close binary system. Interestingly, a TZO nature of WN8 
stars would explain their relatively large stellar radii, high large 
space velocities, the avoidance of OB associations and the (cur- 
rent) dearth of binaries. 
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6. Conclusions 

Summarizing the results of our intense, multi-site monitoring of 
the variable WR star WR 123 in 2003, we note that: 

1 . the star shows pronounced long-term, large-scale, seemingly 
stochastic spectral and photometric variability, broadly in 
line with all previous studies; 

2. the only stable period which is clearly present in the EW 
of prominent spectral lines matches the P = 9.8-hour period 
found in the MOST photometry obtained one year later; 

3. practically all the observed wind volume is involved in the 
10-h periodic variations. Though the variability patterns may 
differ among the lines with different ionization potentials, all 
the lines vary without measurable phase delays. This points 
to the stellar core as sole driver of the periodic variability. 
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